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ABSTRACT 

Giant planet embryos are believed to be spawned by gravitational instability in massive 
extended (R ~ 100 AU) protostellar discs. In a rece nt paper we have shown that dust 
can sediment inside the embryos, as argued earlier bv lBossI (|1998f ) in a slightly different 
model. Here we study numerically the next stage of this process - the formation of a 
solid core. If conditions are conducive to solid core formation, the centre of the gas 
cloud goes through the following sequence of phases: (i) becomes grain (and metal) 
rich; (ii) forms a terrestrial mass solid core via a rapid collapse driven by self-gravity 
of the grains; (iii) starts to accrete a gaseous atmosphere when the solid core reaches 
mass of a few to 10 M®. This sequence of events may build either terrestrial planet 
cores or metal rich giant planets inside the larger gas reservoir of the giant planet 
embryo. In a companion Letter we argue that tidal and irradiation effects from the 
parent star should disrupt the outer metal poor layers of the embryo, releasing nearly 
"ready to use" planets. We propose this as an alternative way to build planets. 

Key words: 



o 
o 



* 



1 INTRODUCTION 

It would be clearly an attractive proposition if planets could 
form in a way similar to stars, e.g., by gravitational con- 
traction followed by a rapid collapse of gas or solid mate- 
rial clumps. A clear route to this is gravitational instability 
(GI) of the parent protoplanetary disc, when the latter frag- 
ments onto gaseous cl umps with masses of a few Jupiter 
mass (Mj = 10" 3 M^I. IBqssI (J1997I ) argued that giant plan- 
ets may form this way. However giant planets in the So- 
lar System contain massive solid cores, and they are also 
much more metal rich than the Sun. Thus a simple gr avita- 
tional collapse of a giant planet embryo could not work. lBossI 
(1 19981 ) suggested that solid cores of the order of a few Earth 
masses (M® ~ 3 x 10 -6 Mq) can be for med if dust grows 
and sediments inside the giant embryos. iBoss et al.l (|2002r ) 
further suggested that even more massive and icy cores can 
be built inside giant proto-planets at 10 to a few tens of AU 
distances. These authors further pointed out that metalic- 
ity of such planets could be raised if gas envelopes of giant 
planet embryos are removed by photo-evaporation due to a 
nearby OB star. Such a model would account for both the 
solid cores and high metalicities of the giant planets. 

However, IWuchterl et al.l (J2000j) argued that giant 
planet embryos may become convective within the first ~ 
100 years of their ev olution, and that co n vection may inhibit 
grain sedimentation. iHelled et al.l (|2008l ) ; iHelled fc Schubertl 
(|2008f ) have recently confirmed this, obtaining much smaller 



solid cores than the model of lBossI (|1998f ) predicts. These au- 
thors also found that their model gas clumps either were too 
hot to begin with or became hot too soon for a significant 
dust growth, instead vaporising the grains. 

Another objection to the gravitatio nal instability m odel 
for planet formation was pointed out bv lRafikovl (| 20051 ) who 
showed analytically that the disc cannot even fragment on 
self-bound gas clumps inside ~ 30 AU at least. The core 
accre t ion (CA) model for planet formatio n (e.g.. IWetherilll 
1 19901 ; iPollack et ail Il996l ; llda fc Linl 120081 ) therefore seems 
more attractive to many. 

However, the extrasolar planets observed at large (R v ~ 
100 A U) separations from their parent stars (jBaraffe et al.l 
|2010j) cannot be explained by the CA model due to ex- 
ceedingly long core formation times. These planets were 
almost certainly formed b y the gravitational disc instabil- 
ity (e.g., see simulat i ons bvlStamatellos fc Whitworthll2008l ; 
llnutsuka et al.ll2009l ; iMachida et alj|20ld ; iRice et aLlboiol ). 

Formation by GI is not absolutely forbidden even for 
the giant planets observed inside a few AU disc. True, the 
planets could not have started to form in situ there be- 
cause initially giant embryos are much less dense than the 
local tidal density (see § 2 below). However, radial m igra- 
tion of massive planets (|Goldreich fc Tremaind Il980l ) due 
to planet-disc torques is now an accepted ingredient of any 
planet formation theory, and provides an attractive expla- 
nation for the observed distribution o f extrasolar plan ets in 
the radial range of ~ 0.1 - 2.5 AU (JArmitagel 120071 ) . and 
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"hot jupiters" even closer in. We do not see why it is not 
possible for the giant planets to form at very large radii by 
GI and then migrate inwards to a few AU if they contract 
quickly enough. In a companion paper (Nayakshin 2010b; 
submitted to MNRAS, referred below as paper III) we ar- 
gue that such a long-range migration is possible in the very 
early phase of formation of the parent star, when the disc 
mass is comparable to that of the star. 

Furthermore, Boss et al, Wuchterl et al, and Helled and 
co-authors all considered embryos with gas densities of or- 
der 10 -8 g cm -3 (perhaps having in mind the tidal density 
constraints for Jupiter at ~ 5 AU). Nayakshin (20 10a, "pa- 
per I") showed the disc fragmentati on constrains IjGammiel 
l200ll : iRice et al.ll2005l ; lRafikovll2005l ) require that the initial 
density of gas clumps formed at ~ 100 AU from the star is 
around 10 -13 g cm -3 . This fluffy initial state of giant em- 
bryos is also much cooler, with typical temperatures of only 
~ 100 K, an order of magnitude lower than in the study by 
Helled an co-authors. 

It is thus not surprising that we found in paper I 
such embryos to be much more pro mising sites o f grai n 
growth and sedim e ntatio n than did iHelled et al.l (J2008I) : 
iHelled fc Schuberj (|2008T ). For reasonable opacity values, 
dust grows and sediments inside isolated giant planet em- 
bryos with masses between a few and 10 — 20 Jupiter masses 
most readily. Since the embryos are approximately isentropic 
at birth, convection did not appear to be important in op- 
posing grain sedimentation. 

For technical reasons, the calculations were stopped in 
paper I when the grain density in the inner region of the gas 
clump reached that of the gas. In this paper we extend the 
calculations similar to those presented in paper I to later 
times. As predicted analytically in paper I, the inner region 
becomes completely grain-dominated by mass and then un- 
dergoes a rapid collapse. The gas component is "protected" 
against the collapse by the pressure gradient, and thus the 
two phases separate at the cloud centre. The grains form a 
high density "solid"' core (it is most likely melted or even 
liquid given the high energy release, but we shall still refer 
to it as a solid core as it may solidify later on when it cools). 

Our main results can be summarised as following: 

• The heat generated by the rapid assembly of the high- 
Z core produces a negative feedback on the core's growth. 
Grains can be melted, driven away convectively or by a 
bulk gas expansion, cutting off further core growth. Thus, 
concerning the final mass of the solid cores , our results 
are le ss optimistic than that of IBose (1199a): iBoss et al] 
J2002J'). but more optimist ic than that of IHelled et al.1 (|200aV 

IHelled fc Schubertl (|2008r ). Continuing the trends from pa- 
per I, the embryos of a few to ten Jupiter masses usually 
yield the heaviest solid cores. More massive giant embryos 
are hotter to begin with and hence reach the grain vapori- 
sation temperature sooner. 

• Gas immediately adjacent to the solid core is strongly 
influenced by the energy release due to the core's assembly. 
Strong convection sets in (thus confirming its importance 
but only after the core formation), frequently supplemented 
by an outward gas expansion. 

• Resolving the region near the solid core is challenging 
numerically, which implies that we probably underestimate 
the degree to which the gas becomes bound to the core. 



Even so, in a simulation where most of the giant embryo 
gets unbound, the solid core is left with a gas atmosphere 
of ~ 0.03 Mjj. This shows that solid cores bred inside giant 
embryos may attract a gaseous atmosphere from within the 
embryo. This atmosphere is metal rich. 

• We then use analytical arguments and a physical anal- 
ogy to the classical results fro m the core accretion model 
for gi ant planet formation (e.g.. lMizunolll980l : |Pollack et al.l 
Il99q ). We suggest that as the solid core mass increases above 
~ 10 — 20 M®, a runaway gas accretion sets in, in which the 
gas envelope collapses on itself. The resulting solid core plus 
a metal-rich self-bound atmosphere should be similar to a 
giant planet. 

• In one of the simulations presented below most of the 
gas envelope gets unbound by the energy release from the 
core. This shows that, at least in the high optical depth 
cases, the heat released by the core's accretion may remove 
the envelope without an outside help. We believe that this 
outer layers blowout process may be even more prominent 
(e.g., likely) when the massive atmosphere collapse occurs, 
as that may release even more energy. 

• We suggest that the outer metal-poor layers of the gi- 
ant embryo could also be removed by physical collisions of 
embryos. 

• We find that the protective giant embryo environment 
strongly diminishes the role of turbulence. We show that 
turbulent gas velocities would have to be as large as 0.3 
times the sound speed to prevent gravitational collapse of 
the "grain sphere" into the solid core. As these high speeds 
are unlikely, we argue that solid cores can be built by a di- 
rect gravitational collapse inside the giant envelopes without 
having to go through the intermediate step of planetesimals. 
The planetesimals are thus not building blocks of planets in 
this picture. They are the material that did not join the 
core. 

• Thus, in contra st to the classic a l planetesimal core as- 
sembly route (e.g., ISafronovl Il969l : IWeidenschillina Il980l : 
IWetheril]||l990n , we would expect a far smaller fraction of 
the solid mass to be locked into intermediate (e.g., a km) 
size bodies. 

• The main building blocks of solid cores in our scenario 
are pebble sized grains (e.g., ~ 10 cm). Grains of this size 
sediment efficiently enough and yet slowly enough as to 
avoid shattering in high speed collisions. 

• Concluding, our key point is that the combination of 
dust sedimentation, solid core accretion and gas envelope 
collapse, on the one hand, with the removal of the outer 
metal-poor envelope of the embryo on the other, may result 
in terrestrial planet embryos, terrestrial planets, or giant 
planets, depending on when the envelope removal occurs. 

In the follow-up paper already mentioned above (paper 
III), we argue that removal of the giant embryo by tidal 
shear from the star helped by heating caused by irradiation 
is un-escapable for giant planets migrating closer in than 
a few to a few ten AU. An exception to this could be a 
more massive giant pl anet M p £ 20Mj that could undergo 
the "second collapse" (|Larsonlll969i : iMasunaga et al.lll99a ) 
before it is disrupted. Such massive planets should have met- 
alicities close to that of their parent stars. 

A hybrid regime of planet formation is possible, in which 
the solid cores are hatched inside the giant planet embryos, 
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delivered into the inner disc and left there when the embryos 
are disrupted. These solid cores may continue their growth 
by accumulation of solid material from the disc (rather than 
the parent embryo) in direct impacts, and then accrete mas- 
sive gaseous envelopes from the surrounding gas disc, result- 
ing in gas giant planets. In this case the beginning phase 
of planet formation is as in the modified GI hypothesis we 
study, but continuation as in the core accretion theory. 

Future detailed calculations and comparison to obser- 
vations are needed to test feasibility of these ideas. In any 
event it seems very clear that the gravitational instability 
model for planet formation has been discounted by most 
authors without a proper critical review of all the evidence. 



It was shown in paper I that grain sedimentation 
timescales can be shorter than the "vaporisation time", t vap , 
which is the time needed for the first core to heat up to the 
grain vaporisation temperature of about 1400 K. The calcu- 
lations were stopped when the density of the dust was equal 
to that of the gas in the inner part of the cloud. Continua- 
tion of the calculations leads to the dust density increasing 
even further, and then a gravitational collapse. This then 
requires a treatment for the solid core forming in the inner- 
most part of the gas cloud. Before we present calculations 
that now include formation of the cores, we first discuss the 
order of magnitude of physical characteristics of the solid 
cores. 



2 PRELIMINARIES 

2.1 Initial properties of self-gravitating clumps 

In paper I we argued that the properties of the giant planet 
embryos formed by gravitational instability in the disc at 
their inception should be similar to that of the first cores 
of same mass. The key properties of first c ores, taken from 
the results of iMasunaga fc Inutsukal (|1999T ) as explained in 
paper I, are briefly reproduced here for reader's convenience. 
The opacity law is taken in the form 
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where kq is the opacity at T = 10 K, and reasonable values 
of a are thought to be between 1 and 2. For simplicity, in 
this paper we only consider the a = 1 case ( H3.5p . The initial 
radius of the first core is given by 



R Cc = 6.0 AU m" 173 ^ 9 " k 4 J 9 



(2) 



where k» = ko/0.01 and T\ — Ti n i t /10K is the initial 
gas temperature, e.g., the temperature of the parent gas 
material (the disc or the molecular cloud). Finally, mi = 
Mf c /10Mj. The mean density of the first core is defined by 
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The mean temperature of the first core is estimated as 



T= - 



1 GM fc /i 



= 146 K m\ /3 T[ 
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3 fcs-Rfc 
The cooling time of the first core is found to be 

icooi = 380 years m/ T^~ kJ . 



(3) 
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(5) 



Note that the cooling time becomes longer with time, as the 
gas clump contracts and becomes optically thicker (paper 

I)- 

2.2 Dust sedimentation in isolated first cores 

The calculations in paper I were done under the assump- 
tion that the first core gains or looses no mass from the 
surroundings. This limit applies if the mass exchange rate 
satisfies (for a = 1) \dM ic /dt\ < 10~ 2 Mj yr" 1 , which does 
depend on parameters of the first core (see paper I). We 
shall continue to assume this in the present paper as well, 
and hence the mass of the embryo is a free parameter of our 
calculations. 



2.3 Maximum masses of solid cores 

The first question to ask about the solid cores built inside 
isolated first cores is how massive they could be. The total 
mass of "high-Z" elements, e.g., those excluding Hydrogen 
and Helium, is about 

M z , t = f g M {c « 60 M e ^ rm , (6) 

where f g = 0.02 is the metal mass fraction (0.02 for gas of 
Solar metalicity). This mass includes ices and organics that 
are probably too volatile to contribute to the solid core, as 
the surrounding gas temperature is typically around 1000 
K. The mass of refractory material (e.g., silicates) that can 
be used to build a massive core is hence expected to be 
about ~ 1/3 of the total heavy elements mass. Therefore, 
the maximum mass of the solid core that can be built inside 
an isolated first core is 
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The minimum core mass (same as the opac ity fragmentation 
limit iLow fc Lvnden-BelJlilffll lRees!ll97cJ ) is a few Jupiter 
masses, and the maximum mass that permits grain sedimen- 
tation is typically a few tens of Mj. Accordingly, the total 
mass of useful core-building material is expected to be in 
the range of ~ 7 to ~ 60 or so Earth's masses. These es- 
timates should be tripled if the embryo is cool enough to 
permit formation of ices. 



2.4 Binding energy of solid cores 

We expect the core of high-Z elements to have a density 
of the order of the terrestrial planet densities, e.g., p c ~ a 
few g cnr 3 . The corresponding radial size of the solid core, 
^?corc = (3M c /47rp c ) , is in the range of 10 8 to a few times 
10 9 cm, which is about 5 orders of magnitude smaller than 
the size of the first core. The binding energy of the solid core 
is 



_ 3 GM corc 

-C/bind.c — _ p 

O Jtcorc 



1.3 x 10 3a erg 



M ffl 



5/3 



-1/3 



(8) 



The binding energy of the first core at formation is in the 
rough range of 10 40 to 10 42 erg (see paper I), and it increases 
in proportion to the virial temperature of the core (dimen- 
sionless T(t) introduced in paper I) as the core contracts. 
Hence a massive enough solid core, M v ~ 10 M®, may re- 
lease an amount of energy comparable to the binding energy 
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of the core itself, strongly affecting the core. Furthermore, 
for a constant density gas sphere, inside the embryo, the 
binding energy scales as M (R)/R oc R , and the thermal 
energy scales as M(R) oc R 3 . Thus building up a massive 
solid core is very likely to heat up the gas at least in the 
inner parts of the first core. 



3 NUMERICAL METHOD 

3.1 Overall approach 

We use the same set up, equations and numerical techniques 
as in paper I, adding to the code several new features to 
describe the evolution of the solid core and its effects on the 
surrounding gas. In brief, spherically symmetric Lagrangian 
coordinates are used. A cell is defined by its fixed gas mass, 
whereas the radial position of the shell changes as the gas 
contracts or expands. The equations are those widely used 
in stellar evolution calculations, except we employ a time- 
dependent approach that does not assume the hydrostatic 
or the energy balances for the gas. We find that this time- 
dependent approach is very important for the problem at 
hand. Both the grain and the gas distributions may evolve 
very rapidly when a heavy-Z core is formed via a rapid (and 
hence quite luminous) accretion of grains. 

The grains are described as a second fluid that is under 
the influence of gravity and gas drag. Grains are allowed to 
slip through the gas, e.g., from one gas radial shell into the 
neighbouring one. For example, grains can move from an 
outer into an inner radial shell when they sediment, or be 
carried in the opposite direction by turbulent or convective 
mixing. 

3.2 Convection and grain convective mixing 

The transport of energy is primarily due to the classical ra- 
diation diffusion flux, although convective flux does become 
important during certain stages of the evolution, e.g., when 
the liquid core's luminosity is high. The mixing length the- 
ory fo r convection is used (e.g., §7 in lKippenhahn fc Weigerj 
ll990T I. We impose an upper limit to the convective flux, F cony 
given by F Cl max = c 3 P, where c s and P are the gas sound 
speed and pressure, respectively. 

In contrast to paper I, the convective motions of the 
gas are allowed to contribute to mixing of grains, imped- 
ing grain sedimentation. Treating this process as a diffu- 
sion process, similar to the turbulent mixing (paper I and 
iFromang fc Papaloizoull2006r ) , we "add" the convective mix- 
ing to turbulent mixing by increasing the turbulent mixing 
coefficient a d introduced in paper I: 



3.3 The inner grain boundary: a heavy elements 
core 



a d = a d + 



-Tcon 



(9) 



This prescription is motivated by the fact that the convective 
mixing length by definition becomes about equal to the gas 
pressure scale-height, e.g., the radius of the first core, R{ c , 
when F con = F Cimax . Note that the maximum of a' d ~ 1. 
Strong convective motions could therefore be quite impor- 
tant in opposing and perhaps even reversing grain sedimen- 
tation, but we find that this fortuitously occurs "too late", 
when a massive solid core has already formed (see [J4] below 
for an example). 



Like previous authors fe.g.. lHelled et alj|200a ). we found it 
very challenging numerically to resolve the radial scales of 
the heavy-Z core for the entire duration of the simulations. 
As there is very little gas mass on these scales, except if and 
when a massive gas atmosphere builds up, many more radial 
shells are needed. These shells require tiny time steps despite 
being not very interesting for most of the calculation. There- 
fore, due to numerical constrains, the length scales that we 
do resolve for most of the simulations are still large com- 
pared with the solid core, of the order of ~ 10 11 cm. Note 
that a snapshot re-simulation can of course be done at a 
higher resolution if one is interested in the small scale struc- 
ture. 

Given that the smallest scales are not properly resolved, 
there is an uncertainty in choosing the boundary conditions 
for grains at R = radius, that is in the very first gas mass 
zone. In paper I we assumed that the grains in the innermost 
zone are suspended by turbulent mixing until the grain den- 
sity in the zone exceeded that of the gas. Presumably the 
inverse drag effects imposed by grains on the gas stifle the 
turbulence once the grains dominate the gas by mass. We 
stopped the simulations at that moment in paper I, delaying 
a study of core accretion to this paper. 

There is no a priory reason why turbulent mixing must 
be effective at R — > 0. If it vanishes in that region, then 
grains may start to sediment into a high-Z core even be- 
fore they overwhe lm the gas by mass in the first zone. 
I Helled et al.l (|2OO80 classified a similar approach as "case 
1", whereas not introducing a solid core until the grains 
dominate the central zone is is analogous to their case 2. 
However, numerical experiments that we did showed that 
the difference in results between these two cases is minor for 
most of the interesting parameter space. Namely, we found 
that before the first zone is dominated by grains, the high-Z 
core grows very slowly in "case 1" anyway. This is due to 
two factors. First of all, the grain density in the central zone 
is small during these "early" times. Secondly, the grains are 
themselves small, yielding a rather small sedimentation rate. 

In the simulations we explored, there is a period of a 
very rapid (but limited) growth of the solid core in both case 
1 and 2. The only significant difference arises for rapidly 
contracting first cores that reach grain vaporisation temper- 
atures before the rapid high-Z core growth phase. In case 1, 
a small high-Z core may be built, whereas it is completely 
absent in case 2. We shall keep this fact in mind, and pro- 
ceed to use only the case 1 in this paper, since it turned 
out to be numerically more stable during the episodes of the 
rapid high-Z core growth. This does not lead to any loss of 
generality of our results. We still cover the appropriate pa- 
rameter space of the problem while avoiding presentation of 
minor differences from very similar runs. 

Therefore, our prescription for the core accretion rate 
is 



dM a 



Md,l{U2 — Ua2)/Rg2 if U 2 — U a 2 > 



(10) 



where M d ,i is the mass of the dust in the first gas zone, R2 
is the outer radius of the first gas zone, and (u 2 — Ua.2) is 
the velocity with which grains arrive into the first zone from 
the second gas zone. For simplicity, we assume that the core 
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has the same material density as the grains, p a . Therefore, 
the outer radius of the core is i? C orc = (ZM cor:l , / Ait pa) 1 ■ 
This radius becomes the new inner boundary for the gas: as 
the core builds up, the inner radius of the first Lagrangian 
mass zone becomes Ri = -Rcore- As the radius of the heavy 
elements core is initially orders of magnitude smaller than 
R2 this does not lead to any significant perturbations for 
the gas. 



3.4 Radiation from the forming core 

Accretion of grains onto the surface of the solid core gener- 
ates accretion luminosity equal to 



^acc — 



dEbin 



GM COIC M, 



corc i «corc 



dt 



lie 



(11) 



The radiative flux boundary condition is modified from the 
zero flux one used in paper I, to read 



Erad(-Rl) 



47T_Rf 



(12) 



3.5 Fiducial parameters 

The outcome of our grain sedimentation calculations depend 
on many parameters, such as: the initial grain size ao; grain 
mass fraction relative to gas, f g ; turbulent mixing coefficient 
ay; maximum velocity « maj[ for grain growth; the opacity 
power law index a; the opacity coefficient no; the mass of 
the first core, Mf c ; and the initial (same as ambient) gas tem- 
perature Tinit. The parameter space is too large to explore 
completely in one paper. 

Fortunately, we find that some of the parameters listed 
above change the end results very little except near regions 
delineating critical behaviour changes. For example, the ini- 
tial grain size is irrelevant in most cases as small grains grow 
very quickly in our prescription due to Brownian motions. 
Grains also evaporate very quickly once gas temperature ex- 



This leaves three parameters that we found to be instru- 
mental in determining the outcome of the calculations, and 
we shall vary these in a relatively broad range: the opac- 
ity coefficient Ko, the mass of the first core, Mf c , and the 
turbulent mixing coefficient, ad- 

3.6 Simulations Table 

The parameters of the simulations selected for presentation 
in this paper, and some of the more important results of 
these simulations are listed in Table 1. The first column is 
the simulation "ID", which starts with a letter M, H or L. 
These stand to indicate the mass of the first core: "H" for 
heavy, M fc = 20Mj; "M" for medium, M fc = 10M, T ; and 
"L" for light, Mf c = 5Mj. The digit next to the first letter 
refers to the opacity coefficient no of the simulation, such 
that Jk satisfies Jk = — log 10 no- For example, ko = 0.01 will 
have Jk = 2 in its ID. Similarly, the symbol a followed by 
a digit J a signifies the magnitude of the turbulent viscosity 
coefficient, ay, so that J a = — log 10 ay. 

The results part of Table 1 starts with two columns 
for two time variables, i cn d and t acC) both in units of 10 3 
years. The former is the duration of the simulation. The 
simulations were run for as long as there were an appreciable 
change in the results or until the simulations stalled due to 
a vanishingly small time step. The latter time entry, t acc , is 
the time of the "runaway solid core accretion phase" . The "- 
" symbol in this column indicates that the simulation never 
entered this phase because grains vaporised earlier than they 
sedimented. The "NA" (not applicable) entry refers to the 
cases when the core accretion is very gradual, so that no 
runaway phase can be defined. This happens when turbulent 
mixing is significant. 

The last three columns in the Table relate to the masses 
of the first cores in units of Earth masses. M corc is the final 
mass of the core, while M c i and M C 2 are its masses at two 
intermediate times, t — 2 x 10 3 and t = 10 4 years, respec- 
tively. These two masses are given to indicate the speed with 
which the cores are assembled. 



case as well. Only when the other parameters of the prob- 
lems are such that the grain growth time, t ST , is about the 
vaporisation time, i vap (see paper I), the exact initial value 
of ao is important. 

We also find that there is a degeneracy in the role of 
some of the parameters. As the interior temperatures of the 
first cores are typically of order of 1000 K, the typical opacity 
of the first core is k ~ ko10 2o \ A larger value of a can be 
compensated for by a smaller value of no to yield a similar 
opacity through the first core, although there may still be 
important differences between these cases near the high-Z 
core or on the outskirts of the first core. 

We shall therefore fix some of the less influential pa- 
rameters at "reasonable" values, motivated by observations 
whenever possible, and also try to avoid model degenera- 
cies. In line with this, we set f g = 0.02, Ti n i t = 10 K, 
a — 1, ao — 1 cm, u max = 1 m s . A logarithmic gas 
mass grid is used, so that the mass of grid cell increases 
with index i as AMj+i = 1.08AM;, where i — 1,2, ...,i max . 
We used imax = 170 cells for most of the runs below, which 
results in the mass of the innermost gas cell of ~ 10 -3 Mgj 
for M fc = 10Mj. 



4 10 Mj FIRST CORES 

Here we present several calculations for first cores with mass 
of WMj. IDs of these runs start with letter "M" in Table 
1. We shall first look at several tests with the same level of 
turbulent mixing parameter, a = 10 -3 , labelled «3 in Table 
1. 



4.1 Medium-high opacity case (Mla3) 

We start with a relatively high opacity case, ko = 0.1. At 
the typical first core's interior's temperature, we have k(T = 
10 3 K) = 10 for this particular choice of kq. This is quite high 
but may well be reasonable at least in the innermost parts 
of the first core when dust sediments there. 

Figure Q] shows the evolution of several key quantities 
as a function of time for this calculation. The solid and dot- 
ted curves in the upper panel of the figure show the grain 
size, a, and central gas density, p, in units of 10 -9 g cm -3 , 
respectively. The dashed curve is the central gas tempera- 
ture in units of 100 K. The solid curve in the lower panel 
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Figure 1. Results of simulation Mla3 (see il4.ll and Table 1). 
Upper panel: Grain size a in cm (black solid curve), dust-averaged 
temperature T m (dashed green) in units of 10 3 K, central gas 
density in units of 10 — 9 g cm -3 (blue dotted). Lower panel: mass 
of the high-Z elements core, in Earth's masses (black solid), core's 
accretion rate in 10 — 2 M^g yr -1 (red dashed), and core's accretion 
luminosity in 10 30 erg s" 1 (black dotted). 



shows the mass of the high-Z core in units of Earth's mass, 
whereas the dotted and the dashed curves show the core's 
luminosity (in 10 30 erg s _1 ) and accretion rate, in 10 -2 M ffl 

yr" 1 . 

The first ~ 5 x 10 3 yrs of the simulation are familiar 
from paper I. The grains must grow to ~ 10 cm size before a 
significant dust sedimentation takes place. The evolution af- 
ter that is marked by two stages: the "runaway" stage during 
which the high-Z core grows very rapidly, almost vertically 
in Figure [T] and the saturated growth stage when the accre- 
tion rate, core's luminosity and grain size all drop with time. 
This second stage is a testament to importance of feedback 
effects during the core's growth. This can be deduced by not- 
ing the abrupt central temperature increase at the end of the 
"runaway phase", when both the high-Z core mass and lu- 
minosity increase sharply, and the corresponding central gas 
density decrease. The temperature increase disrupts the hy- 
drostatic balance in the central regions, causing an outward 
gas expansion. 

These effects are seen in Figure [5] that presents the ra- 
dial structure of the first core and the grain density distribu- 
tion during and just after the end of the runaway core growth 
phase. The panels show the gas and the grain densities, the 
gas temperature, the radiative (black curves) and the con- 
vective luminosities (blue curves), and the gas (black) and 



the grain (red) radial velocities. The convective luminosity 
is defined by L con = 47ri? 2 F CO n. 

The first snapshot in figure [2] shown with the solid 
curves, corresponds to the time when the grain density in the 
inner region is equal to that of the gas, which corresponds 
to when we stopped the simulations in paper I. The mass 
of the high-Z core is only ~3x 10~ 4 Mj at that point. As 
the grain density overtakes the gas density, the gas drag and 
turbulent mixing resistance to grain sedimentation rapidly 
weaken. This is seen in the inner part of the radial veloc- 
ity curves sequence. As the inner part of the "grain sphere" 
contracts due to its own self-gravity, the accretion rate onto 
the high-Z core increases by orders of magnitude. The lumi- 
nosity of the solid core shoots up to about a tenth of that 
produced by the whole first core (see the dashed and the 
dot-dashed curves in the upper right panel of the figure). 

The heat produced by the build up of the high-Z core 
strongly affects the inner region of the cloud, heating up the 
gas there. The increased gas pressure initiates a rarefaction 
wave that moves gas outward, as can be seen in both the 
density and the velocity curves. It is likely that the grain 
sedimentation would continue at this stage if it were not 
for the gas entropy profile becoming strongly unstable to 
convection. The convective flux builds up to its maximum 
F Cimax value allowed (blue curves in the upper right panel) in 
the inner part of the cloud. Strong convective mixing drives 
some grains in the inner region outwards (see the red veloc- 
ity curves for grains). Note a very good correlation in the 
positions of the peak of gain velocity curves in the inner 
part of the embryo with the maximum in the convection lu- 
minosity. This confirms that the outward motion of grains 
is driven by convection in this region. 

The second, also quite important way in which the high- 
Z core growth self-regulates is vaporisation of grains. As the 
inner region of the computational domain heats up, grain 
growth and evaporation start to compete with each other. 
In fact the grain size decreases with time (upper panel in 
fig. [TJ , although not particularly rapidly in this simulation. 

The existence of the two ways in which the high-Z cores 
can affect their surroundings and hence limit the rate of their 
own growth suggests that this feedback loop is quite robust 
qualitatively, despite quantitatively depending on details of 
the gas-grain interaction prescriptions and the parameter 
values. 

These feedback effects preclude incorporation of all of 
the available condensible material in the high-Z core. Indeed, 
as much as ~ 20 Mg of such material is present in the first 
core, but only a quarter of it gets locked into the high-Z 
core. Clearly this result must be sensitive to the opacity of 
the first core, as lower opacity would allow more heat to 
escape the gas cloud, perhaps reducing the feedback effects. 

4.2 Medium opacity case 

We now present calculation M2a3 which is identical to Mla3 
except the opacity is reduced to kq = 0.01. The result- 
ing evolution of the first core and the grain distribution 
is shown in Figure [3] Confirming the earlier made claim 
that embryo opacity is a key parameter of the problem, the 
resulting change in the dust sedimentation outcome is sig- 
nificant. As the cooling time of the first core decreases as 
Ko decreases, the cloud contracts more rapidly. Higher gas 
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Figure 2. Radial structure of the first core during and just after the "runaway growth phase" in the simulation Mlo3 (see Fig. [TJ. 
Densities and velocities of gas and grains are shown with the black and red curves, respectively. The different style lines are for times 
t = 5386, 5417, 5446 and 5469 years, for the solid, dotted, dashed and dash-dotted curves, respectively. The blue curves in the right upper 
panel show the convective luminosity which becomes relatively large after the high-Z core creation. 



density promotes faster grain growth. Therefore, grains in- 
crease in size and sediment sooner in M2a3 than they did 
in simulation Mla3. The runaway accretion phase begins 
earlier, at around 2 x 10 3 years. As the core accretion rate 
is higher, the feedback effects of the growing high-Z core are 
more pronounced, and hence the rapid growth of the core 
saturates at a lower core mass, and a lower luminosity (cf. 
the lower panels of figures [T] and |3J. As in simulation Mla3, 
the gas temperature near the core is above vaporisation tem- 
perature for the grains, causing grain melting and hence a 
plunge in the value of a at the peak of the runaway growth 
phase. In comparison to Mla3, reduced opacity enables the 
gas near the core to transfer the heat outward more rapidly. 
This allows a further density and temperature increase near 
the core, evaporating the grains further. A complete grain 
evaporation occurs at about 3,500 years in this calculation. 
Grain accretion onto the high-Z core then ceases, allowing 
the gas near the core to cool slightly. The contraction and 



the respective heating of the gas near the core continues. 
This stops any further core re- growth. 

The final mass of the high-Z core in this simulation is 
M core = 0.9996 M e - The fact that it so close to the Earth's 
mass is of course a complete coincidence. Degrading the nu- 
merical resolution from i ma x — 170 to i ma x — 140, and 
repeating the simulation M2a3 leads to qualitatively very 
similar results but the final core mass is M COIC = 1.24 M®. 
This implies that higher numerical resolution calculations 
would yield somewhat lower masses for the final M core , but 
probably not by more than a factor of 2. This is acceptable 
to us given large uncertainties in parameter values. The de- 
crease in the final core mass with an increasing resolution is 
probably due to feedback effects becoming even more pro- 
nounced as regions closer and closer to the core are resolved. 
We shall investigate these issues further in the future. 

Figure U shows the radial structure of the first cores at 
and immediately after the runaway accretion phase for the 
simulation M2a3. There are many features similar to those 
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Figure 3. Same as figure ITl but for simulation M2ct3. Note the 
anti-correlation between the size of grains and the central tem- 
perature after the high-Z core is created. This is due to a feedback 
loop: grains power core's accretion luminosity that heats up the 
central region of the first core, but the grains can then be melted 
by if the gas becomes too hot. 



seen earlier in figure [2] Note a slightly lower core luminos- 
ity, which is in line with the lower final core mass. It is also 
interesting that the general grain sedimentation trend is re- 
versed at "large" R, R ~ 2 x 10 12 cm, as grains there start 
to move outwards (see the lower right panel in figure |4]). 
This is caused by a rapid decrease (by a factor of two) in 
the grain size at the runaway accretion phase. The reduced 
size of grains implies a smaller sedimentation velocity, and 
sedimentation is apparently overpowered by turbulent mix- 
ing at these radii. The effect however does not last for long 
as grains re-grow and sedimentation continues. 



4.3 Low opacity case 

Simulation M3a3 probes the case of a yet smaller opacity, 
K — 0.001. The cooling time of the first core in this case is 
very short. Therefore, the results of this test are "uninter- 
esting" - grains grow to a m 1.5 cm only to be evaporated at 
about t m 200 years. We present no figures for this calcula- 
tion. Shortly after the grains evaporate, the gas temperature 
approaches 2000 K, and the first core shoul d undergo the 
second collapse (jMasunaga fc Inutsukall2000l ) and probably 
form a giant planet mass object without a solid core. 



Simulation M0a3 is the highest opacity case we consider 
here, with ko = 1. As the cooling time is very long for this 
simulation, the first core does not contract very much be- 
fore the grains sediment into the high-Z core. The central 
density of the first core is much smaller than in the pre- 
viously investigated cases, and for this reason it is plotted 
in units of 10~ 12 g cm -3 in Figure [S] As can be seen from 
the figure, the simulation M0a3 is very different from the 
lower opacity cases studied so far. The central temperature 
of the first core in M0ev3 never becomes high enough to evap- 
orate the grains, except at the very end of the simulation. 
This is why the grain size (solid curve in the upper panel 
of figure [SJ) and the accretion rate onto the solid core (red 
dashed curve in the lower panel of the figure) keep increas- 
ing through and after the runaway accretion phase, whereas 
in Mla3 and M2a3 these quantities peaked during the run- 
away phase and subsequently dropped with time. 

Due to lower gas temperatures, the melting type of feed- 
back does not take place in M0a3 before t~3x 10 4 yrs. 
Grains continue to sediment and build up the central solid 
core at an increasing rate. The heat released by the core is 
trapped near the core due to a higher opacity of the embryo. 
The only way for the gas in the inner part of the first core 
to absorb the energy released by the high-Z core is to put it 
into a kinetic form. This drives a strong outward expansion. 

Figure [6] shows the radial structure of the first core near 
the beginning of the high-Z core runaway accretion phase 
(solid curves) at t = 8.5 x 10 3 yrs, and at a later time, 
t — 2.44 x 10 4 yrs (dotted curves). Note an immense change 
in the luminosity curves (upper right panel of the figure). 
While initially heat production was all due to the compres- 
sional heating of the gas, at later times the energy produc- 
tion is dominated by the grain accretion onto the high-Z 
core. The innermost region of the gas becomes unstable to 
convection, and the convection flux then increases to the 
maximum sonic value that we impose ( H3.2p . At the same 
time, the convection flux is still very small compared with 
the radiative flux. 

The increased luminosity of the high-Z core also forces 
the inner regions of the gas to become much hotter, although 
not as hot as to evaporate the grains (see the lower left panel 
of figure [6}. The rarefaction in the inner ~ 10 12 cm of the 
first core is significant, e.g., gas density decreases by a factor 
of a few. Note that the motion of the gas and the increased 
convective grain mixing does cause some outward radial mo- 
tion of the grains, although the innermost grain-dominated 
region becomes depleted mainly due to grain accretion onto 
the high-Z core (see the red curves in the upper left panel 
of the figure). 

Further evolution of the first core is marked by a contin- 
uing energy release from the solid core, which continues to 
inflate the hot bubble within the first core until it actually 
breaks out - unbinds most of the gas in the outer parts of 
the first core. This can be seen in Figure [7] which presents 
the radial structure at time t — 3.16 x 10 4 yrs. Note the 
very high temperature and density at small radii and, on 
the contrary, very low densities and a fixed temperature of 
20 K at large radii. The latter is the minimum temperature 
imposed at the outer boundary for this simulation (and for 
most others; normally this plays no role at all as the tern- 
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Figure 4. Same as figurc[2]but for simulation M2«3 (see also figure [3). The snapshot times are t = 1880, 1892, 1902 and 1912 years. 



perature of the outermost zone is usually higher than 20 K). 
The radial structure of the gaseous part of the first core in 
figure [7] can be divided into three parts : (i) the dense and 
hot envelope tightly bound to the high-Z core; (ii) the inter- 
mediate region at 10 12 £ R ^ 10 13 ' 5 cm; and (iii) the outer 
isothermal expanding shell. The last part of the fist core is 
the most massive one; it is completely unbound and is lost 
to the environment. 



On the other hand, the region near the high-Z core can 
be called its atmosphere as it is apparently bound to the 
core. It is likely that the removal of the outer massive shell 
of gas actually promoted formation of the atmosphere by 
letting the heat escape the inner region more easily. The 
atmosphere is not massive at all: the upper left panel of 
figure [7] shows that it is around ~ 0.03 M®, compared with 
the mass of the core, M corc ~ 20 M®. Nevertheless, it is 
important as it indicates that the solid cores formed inside 
the first cores can actually build up gaseous atmospheres if 
right conditions exist. 



4.5 Nonlinear effects of turbulent mixing 

We found that the models can be surprisingly and non- 
monotonically sensitive to the turbulent mixing parameter 
ad- Simulations M2q2 and M2a4 are identical to M2a3 ex- 
cept that a c i = 10~ 2 and 10~ 4 , respectively. The evolution 
of the first core properties of simulation M2q?2 and M2q4 
is plotted in Figure [Sj the left and the right panels, respec- 
tively. 

Concentrating first on M2a2, we observe that its out- 
come is drastically different from that seen in figure [3] The 
larger turbulent mixing delays grain sedimentation by only 
about 50% (see i acc in Table 1), but the gas in the first core is 
hotter by this time. Importantly, the gas temperature is now 
closer to the vaporisation temperature. As is clear from the 
upper panel of figure [8^,, accretion luminosity of high-Z core 
is only ~ 10 27 erg s — 1 , butthis is sufficient to push the gas 
near the core over the vaporisation temperature threshold. 
The grains evaporate rapidly soon after this, and a further 
growth of the high-Z core turns out impossible. The final 
mass of the core is M coro = 0.08 M®. Note also that the 
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Figure 6. The radial structure of the first core near the end of simulation M0a3 (see Fig- [5j. The solid and the dotted curves correspond 
to times t = 8.5 X 10 3 and t = 2.44 X 10 4 yrs, respectively. The red dotted curve in the top right panel shows the convective flux, 
saturated at the sonic value in the innermost part of the first core. 



accretion rate of grains on to the high-Z core is smaller in 
M2a2. 

The simulation M2a4 has a much lower level of turbu- 
lent mixing. Recalling the trend from M2a2 to M2a3, e.g., 
the increasing mass of the solid core with the decreasing 
turbulence parameter, one could expect the solid core to be 
more massive in M2q4 than that in M2a3 (1 M®). However, 
the result is actually M coro ~ 0.7 Mq. The right panel of 
Figure [H] shows that the time evolution of the central part 
of the cloud in this simulation is quite complex. The max- 
imum core accretion luminosity in this simulation actually 
exceeds that in M2a3 (compare the sharp red spikes in the 
luminosity curves for these two simulations). The energy re- 
leased by the core then melts the grains rapidly. While in 
M2a3 grain re-growth occured almost immediately after the 
accretion luminosity drops, in M2a4 the gas near the core 
is a little hotter and hence the grains remain small for sev- 
eral hundred years. After the excess heat finally leaves the 
central region, grain growth restarts. This apparently leads 
to a second major core accretion episode at around 2200 



years. The intricate feedback loop steps in again, melting 
the grains somewhat, but this time to the size of about 1 
cm rather than ~ 0.1 cm. The grains are eventually vapor- 
ised at time t = 3300 years, very similarly to the simulation 
M2q3. The core mass is then frozen at the final value of 
0.7M ffi . 

Comparing the three simulations, it appears that the 
role of the turbulent mixing is to smooth out the evolution 
of the grain size and the core accretion rate. Indeed, we note 
that turbulent grain mixing, in prescription we are using, 
may have both a negative and a positive effect on the solid 
core growth. Prior to the formation of a relatively massive 
core, turbulent mixing works only in one direction, that is 
opposing grain sedimentation. However, once a sufficiently 
massive core forms, gravity in the inner regions become dom- 
inated by the solid core itself rather than by the gas. The 
sedimentation velocity of the grains then increases above the 
"pre-core" theoretical value (paper I). The inner few zones 
of the core are then quickly emptied of the grain material by 
accretion onto the solid core. This has the effect of creating a 
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Figure 7. The radial structure of the first core near the end of simulation M0«3, at time t = 3.16 X 10 4 yrs. The upper right panel 
shows the enclosed gas mass in Earth's masses. 



positive rather than the usual negative gradient in the ratio 
of Pd/p, e.g., a relative deficit of grains in the centre of the 
first core. Turbulent diffusion mixing then reverses its effects 
and actually increases the grain accretion rate near the core 
while of course continuing to impede grain sedimentation at 
larger radii. 



5 DEPENDENCE OF RESULTS ON THE 
MASS OF THE FIRST CORES 

Having studied grain sedimentation and high-Z core building 
for Mf c = 10Mj in some detail, we shall now review the 
results of the less and more massive runs shown in Table 1. 



5.1 "Heavy" first cores 

The most massive first core case studied in the paper is 
Mf c = 20Mj, shown in Table 1 with IDs starting with letter 
"H". Only one value of turbulent parameter is considered 
for these simulations, ay = 0.001. It is immediately obvious 



that these most massive cores are the least promising sites of 
grain sedimentation, which is also in accord with the results 
of paper I. Only the high opacity cases, kq = 0.1 and 1, 
resulted in significant high-Z cores. The physical reason for 
this is their already hot initial state. These cores rapidly 
evaporate their grains, if allowed to cool at a reasonably 
rapid rate. 

The inefficiency of the process of dust sedimentation for 
these most massive cores stands out clearly. The most mas- 
sive solid core assembled, M corc = 12.2 M®, represents only 
about a quarter of the condensible high-Z mass available 
inside the gas clump for Mf c = 20Mj. 

5.2 "Low mass" first cores 

Finally we come to the lightest first core case we studied 
in this paper, Mf c = 5Mj. Concentrating on the first four 
entries starting with letter "L" in IDs in Table 1, one notices 
that these runs all produce high-Z cores with masses above 
1M®. The efficiency with which the grains are locked into 
the central core reaches 100% in L0a3 and Lla3. 
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Figure 8. Same as figurc[T]but for simulation M2o2 (Left panel) and M2a4 (Right panel). Note that in both simulations the end mass 
of the solid core is lower than that in A2a3, but for different reasons. See text f £|4.51 for a further explanation. 



Looking at further entries in the Table, we find the 
familiar trends with the turbulent mixing parameter ad' 
higher values of it usually result in smaller values for M corc . 
The turbulence needs to be rather extreme to prevent the 
formation of high-Z cores completely: a very high turbulence 
run, ay = 0.1 (LOal), resulted in a core of 0.075 M^ mass, 
albeit after a very long (~ 10 years) time. 

The physical reason for the resilience of dust sedimen- 
tation process in lower mass cores is that the initial temper- 
ature in such is low, e.g, ^ 100 K (see the lower panel of 
Figure 1 in paper I). Both cooling and vaporisation times 
in the low mass first cores are long, and hence dust has a 
plenty of time to grow and sediment, even in the presence 
of turbulence. 



6 LONG TERM EVOLUTION 

6.1 Expectations: Giant planets with solid cores? 

The solid core - gas envelope structure we study here is 
essentially the same as that studied by many authors in 
the context o f the core accretion model for giant planet for- 
mation (e.g.. IPollack et al.lll996T ). The difference is mainly 
in the environmental factors by which we mean boundary 
conditions "far" from the solid core. Furthermore, the solid 
cores that we study here also grow by accretion of solids. 
We should thus expect that at least some of the classical 



core accretion model results may pertain to the problem we 
study. 

One such result concerns the build up of a massive 
gaseous atmosphere around the solid core. For low solid 
core masses, the atmosphere is in a hydrostatic equilibrium, 
in which the gravitational force due to the solid core and 
the gas own weight i s bala nced by the pressure gradient. 
As shown by iMizund l| 19801 ) , when the solid core mass in- 
creases and exceeds a critical value, Af cr it, a hydrodynamical 
collapse driven by self-gravity of the gas takes place. The 
gaseous envelope contracts to much higher densities. This 
marks transition from a hydrostatic balanced atmosphere 
to a gas accretion phase. The latter goes on as long as there 
is gaseous fuel supplied from outside. Adding accretion of 
solids from the protostellar disc one ar rives to a very crud e 
sketch of the core accretion model (e.g.. |Pollack et al.lll996T ). 

There is a scope for a conceptual confusion here. The gi- 
ant planet embryo is itself a massive gaseous self-gravitating 
"atmosphere" around the solid core. But the bulk of this gas 
is not bound to the solid core at all as the mass of the latter 
is rather small. Initially, the gas pressure gradient is nearly 
zero in the central part of the gas cloud since weight of the 
central region is small (proportionally to ~ i? 3 ). It is the 
weight of the layers above that keep the innermost gas from 
expanding outward. Appearance of a solid core in this con- 
stant pressure, constant density, environment, at first has 
little effect as the gas pressure is too high. However, as the 
mass of the solid core increases, there is a critical point as 
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and density values near the end of the simulation. This marks 
building of a dense and hot atmosphere around the solid core (cf. 
figure [Sj . 



in the calculation of lMizund |l980J) at which the gas around 
the solid core starts to collapse on it and forms an atmo- 
sphere around it. This "atmosphere within and atmosphere" 
is gravitationally bound to the solid core and itself. It is 
much denser and becomes distinct from the giant embryo. 
In paper III we argue that when the embryo is disrupted 
by tidal forces of the parent star, the inner solid-core bound 
atmosphere may survive. The result is probably not unlike a 
giant planet built in the core accretion model, although de- 
tailed calculations beyond the scope of our paper are needed 
to characterise its chemical composition and other proper- 
ties. Nevertheless, one robust prediction is a metal rich com- 
position of a giant planet built inside a more massive giant 
embryo gas clump: the abundance of dust in the central re- 
gions of the embryo is enhanced by dust sedimentation. 

6.2 Numerical issues 

The most straightforward way to explore these ideas would 
be to study them numerically. In the simulations listed in 
Table 1, there is only one example of a bound gas atmo- 
sphere, the simulation M0q3, and even then the atmosphere 
weighs only 0.03 M®, which is about 0.001 of the mass of 
the high-Z core for that simulation. 

We suspect that at least the more massive of the solid 
cores built up in our simulations should be surrounded by 



massive gaseous atmospheres at later times. Computational 
limitations of our method precluded us from studying this 
interesting question in detail. 

As explained in H3.ll we need a hydrodynamics rather 
than hydrostatics based approach to follow the evolution of 
the system during the phases of a rapid gravitational col- 
lapse or the bounces driven by too much accretion onto the 
solid core. Unfortunately this hydrodynamical approach also 
makes it difficult to evolve the system for long periods of 
time after the formation of the solid core if the core starts 
to build up a massive gaseous atmosphere. As an atmosphere 
starts to build up around the solid core, time steps become 
very short, and the code basically stalls. 

Being unable to follow the evolution of the system into 
the appropriate regime numerically, we resort to analyti- 
cal arguments and analogies to the classical core accretion 
model results. 



6.3 Minimum critical (cold) solid core mass 

The temperature profile of the gas is the main uncertainty 
in our calculation here. The lower the gas temperature, the 
smaller is the gas pressure gradient for a fixed density profile. 
Therefore, by choosing appropriately a minimum possible 
temperature that the gas near the solid core may have we 
should be able to obtain the minimum required solid core 
mass for atmosphere collapse. 

A simple way of doing this is to consider a giant embryo 
unperturbed by the energy release from the solid core, set- 
ting Lace = 0. This is satisfied by those of our models that 
had grains evaporated after the core assembly (e.g., see the 
right panel of figure [8]). Further, somewhat unrealistically, 
we shall assume that the gas remains isothermal even in the 
presence of the solid core's potential. The result is neverthe- 
less interesting as a lower limit to the mass of the solid core 
at which a massive gaseous atmosphere can be built. 

With these ideas in mi nd, our tre a tment here echos that 
presented in section 6.1 of lArmitaeel (120101 ) (note that the 
derivation was first made by Sasaki 1989|). The hydrostatic 
balance equation for isothermal gas surrounding the solid 
core is 



dp 
dR 



GM C , 



GM e (R) 



(13) 



where 

l-R 



R 2 R? 

c a is the isothermal sound speed, M e (R) — 
J 4-irr 2 p(r)dr is the enclosed gaseous envelope mass in- 

side radius R, and i? C orc is the radius of the solid core. If a 
hydrostatic balance solution exists, then the density given 
by equation [13] should smoothly join the unperturbed first 
core density at radii R where the gravity of the solid core is 
negligible. 

We first consider the case when the gaseous envelope is 
much less massive than the core, and thus we neglect the 
last term in equation 1131 The solution is then trivial: 

uilicorc 



p(R) - pi C exp 



Rri 



for R> R c 



(14) 



Note that p(R) — > pt c at R — >• oo, as needed. 

The factor in the exponential can be very large near 
the solid core. In that case it also drops off very quickly 
with the distance from the first core. Hence, the mass of the 
atmosphere around the solid core is approximately 
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(15) 



This treatment is valid as long as Af atm <C M core , since 
in the opposite case the weight of the atmosphere cannot be 
neglected in equation ll3l as we have done above. Self-gravity 
of the gas around the solid core initiates a "collapse" of the 
gaseous atmosphere. The critical mass of the solid core when 
the collapse happens is found by setting M atm = M core . We 
call the corresponding core's mass a "critical cold" mass: 



GM a 



lie 



ln[^ 

Pic 



20A 



(16) 



where A = ln(p c /pf c )/20 « 1 for pi c ~ 10 -11 g cm~' A . As 
the first core density enters logarithmically, the result is 
weakly dependent on the value of A. Eliminating -Rcorc = 
(3M corc /47rp c ) 1/3 , we have 
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This expression is similar to that for Jean's mass, but with 
the density of the solid core rather than the gas density 
inserted, and the logarithmic factor A. Numerically, 

3/2 /„ „ r \ 3/2 
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This estimate is robust in the sense that the temper- 
ature of the unperturbed first cores does not vary all that 
much. As found in paper I, by the time grains sediment to 
the centre the temperature of the gas is typically between a 
few hundred K and 1400 K. The upper value for the tem- 
perature range is set by vaporisation of grains, whereas the 
lower one is simply due to a rapid cooling of first cores: as 
we found in paper I, the first cores cool initially very rapidly 
since they are spatially extended and only mildly optically 
thick. 

A correction downward to the minimum critical core 
mass is needed in the case when the density of the dust 
around the solid core is comparable or higher than that of 
the gas. One should then adjust the mean molecular weight 
in equation [TS] to reflect the metal-rich composition of ma- 
terial in the inner region of the first core. 

6.4 Maximum critical (radiative) core mass 

The isothermal envelope assumption that we used above al- 
lowed us to explore one extreme of the problem at hand, e.g., 
when the envelope temperature is completely unaffected by 
the presence of the solid core. An opposite extreme assump- 
tion is that the immediate envelope around the solid core 
is much hotter and denser than the giant embryo in which 
it is embedded, so that the presence of the latter is a mi- 
nor detail. This may nearly be the case if the luminosity of 
the solid core is very high, e.g., comparable with the total 
cooling luminosity of the giant embryo. This set of assump- 
tions constitutes the "r adiative zero" so lution to the gaseous 
atmosphere structure (|Stevensonll 19821 ) . 

The radiative zero solution also makes a prediction for 
the critical solid core mass (see below), which we shall call 
M ra d- In contrast to the "critical cold" mass derived above, 
-Mrad is likely to overestimate the critical core mass. The 
radiative zero solution assumes zero pressure and density 



boundary conditions at infinity. In contrast, our envelope 
is surrounded by material of a relatively high density and 
pressure. At large distances from the solid core, gas density 
and temperature should join the values of these quantities in 
the first core, rather than go to zero. The additional pressure 
and weight of the outer cooler envelope of the giant embryo 
is certain to lower the critical mass of the solid core required 
to initiate the collapse. 

The central quantity in the radiate zero solution is the 
luminosity of the solid core, which is derived from accre- 
tion of planetesimals. We thus need to estimate the typical 
accretion rate onto the core. The sedimentation time scale 
of grain material, t aB d has been estimated in paper I, and 
is usually between 10 3 to a few x 10 4 years. Assuming that 
Mz, c ~ 20 M®, and that all this solid core-building material 
falls onto the core within 10 4 years, we find the solid core 
accretion rate, M cor e, is 
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2 x 10 M® yr" 



(19) 



Ikoma et al. 2000 updated and expanded the results of 
IStevensonl (|l982l ) for a broad range of opacities in the enve- 
lope (which were assumed temperature independent in con- 
trast to our power-law dependence). In particular, they pro- 
vided a power-law fitting formula to the critical core mass. 
Using that formula for a fiducial opacity coefficient value of 
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where q' w s' ~ 1/4. 



6.5 Summary on critical core mass and further 
complications 

We have argued above that there should be a critical solid 
core mass, Merit, above which the initially tenuous gas at- 
mosphere around the core is unstable to self-gravity and 
should collapse. The collapse should build a massive gas 
atmosphere atop the solid core. The composition of the at- 
mosphere should be metal rich as the interior of the giant 
embryo is metal rich due to dust sedimentation. This object 
inside the giant embryo may have properties si milar to gi- 
ant g as and giant icy planets (as pointed out bv lBoss et al.l 
120021) . 

As far as the value of M cr it, we suggested that it is 
limited from the bottom by M co id and from the top by M ra d- 
For typical parameters the range here is from a few to few 
tens of Earth masses. 

However, for the most massive embryos, say Mf c £ 200 
Jupiter masses, it is possible to get the gas to collapse grav- 
itationally and form a compact proto-planet without a solid 
core at all (see below), especially if dust opacity can be made 
very small. In this case the embryo may reach vaporisation 
temperature before dust sediments. 

Our numerical code assumes that hydrogen is molecu- 
lar, which is a very reasonable assumption for the first cores 
before they undergo the second collapse (JMasunaga et al.l 
1 19981 : iMasunaga fc Inutsukall2000r ) . This assumption may be 
violated when a dense and hot gas atmosphere builds up 
around the core, or when the whole first core contracts suf- 
ficiently to surpass T ~ 2000 K. As with collapse of first 
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cores, hydrogen molecules disassociation and then hydrogen 
atom ionisation energy losses at even higher temperature 
may help gravitational collapse of the atmosphere. 

Therefore, if the first core is truly isolated and is al- 
lowed to contract in peace (rather than being tidally or ir- 
radiationally disrupted after a finite amount of time, pa- 
per III), then the whole first core can eventually collapse 
onto the solid core. This "second" wave of collapse is com- 
pletely analogous to that of the first cores in star formation 
(|Larsonlll969l ). and is not at all initiated by the solid core. 
Hence truly isolated giant planet embryos could presumably 
contract into giant planets with arbitrary small solid cores, 
including none. 



7 ON ASSEMBLY PROCESS OF 
TERRESTRIAL PLANET CORES 

In ij6] we proposed that giant planets with solid cores may 
form inside the giant embryos when the solid cores grow 
more massive than the critical mass, M cr . We further argued 
that this part of the planet formation picture may be quite 
similar to that worked out in the classical core accretion 
schemes, albeit for different accretion rates of solid material 
and boundary conditions "far" from the solid core. 

In contrast, the earliest stages of the solid core assembly 
could not be more different. We believe it is far "easier" to 
form massive solid cores inside the giant embryos at large 
distances from the star than it is in situ at R ~ f ew AU in 
the "standard" botto m-up build-up scenario (e.g.,[Safronov| 
1 19691 ; IWetheril]||l990l ) . The numerical treatment of dust sed- 
imentation inside the first cores (= giant planet embryos) 
that we employed in this paper may have obscured this im- 
portant point, which we shall discuss below in detail. 



7.1 The standard bottom-up scenario: the key 
role of planetesimals 

The c ore accret i on model - (e.g. , Safronovl Il969l ; iMizund 
Il98d : IWetherilll Il99d: llda fc Linl 120081 ) . also reviewed in 
chapters 4-6 in lArmitagej ( 20101 ) - stipulates that planet for- 
mation starts on smallest scales, with dust initially growing 
inside the protoplanetary disc into larger grains. This very 
fir st pa r t of t he picture is physically similar to the model 
of [Boss] (|1998r ) and what we assumed here, with dust grow - 
ing by gentle hit-and-stick collisions (|Blum fe Wurmll200SI ). 
There is then a poorly understood step of building planetes- 
imals, which are solid bodies of ^ km size. Once planetes- 
imals are assembled, "planetary embryos" - 10 3 km sized 
bodies are built by collisions of planetesimals, and finally 
Earth mass terre strial cores are accreted by giant impacts 
of such embryos (|Safronovll 19691 ). 

Due to well known problems with planetesimal assem- 
bly authors used to urge a leap of faith, arguing that if 
planetesimals had not been assembled somehow, we would 
not be here, a nd that "Nature knew how to do it" (e.g., 
IWetherilll |l_990j) . This is clearly unsatisfactory. Most solid 
matter outside the Sun in the Solar System is inside the 
cores of the giant planets and also in the terrestrial planets. 
Thus the total mass of planetesimals must have been at least 
equal to the total mass of the solids in the Min imum Mass 
Solar Nebula (MMSN, IWeidenschillingl 1 1977bl ) . and hence 



most of the solids in the disc (= "nebula" ) must have been 
reprocessed through planetesimals. 

Recent conceptual breakthroughs suggest that solids 
may be concentrated into larger structures by instabilities 
and turbulence in th e disc (e.g., Youdin fc Goodman! 12003 ; 
IJohansen et alj|2007l ; ICuzzi et al.ll2008l ). although it is too 
early to say that the planetesimal step assembly is under- 
stood. Studies of collisional evolution of asteroid-like bodies 
suggest that the initial size of planetesimals must be larger 
than the classical ~ few k m size objects, and be i n the 100 
to 1000 km size range fe.g.. lMorbidelli et al.ll2009l ). thus in- 
creasing the requirements for the planetesimal buildup mod- 
els. 

In any event, we shall now review the well known rea- 
sons why terrestrial planet cores cannot be built in the CA 
model directly, and argue that these difficulties are avoided 
if massive solid cores are made inside giant planet embryos. 



7.2 Kepler shear: limiting the mass of fragments 

The most direct route to "large" solid objects in the CA 
model is the gravitational instability of a dense dust layer 
(|Safronov|[l969l ; iGoldreich fc Wardlll973T ). It works best if 
there is no turbulence in the gas, as turbulence provides 
support against gravitational collapse. We hence first as- 
sume that turbulence is absent. Inside a protostellar disc, 
Kepler shear limits the maximum mass of the clump of solid 
material that c an collapse into one body to (e.g., §4.6.2 in 
lArmitagdlJHOT l 
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where E so i is the column density of solid material in the disc, 
and Q is the Keplerian ang ular frequency at that location. 
At 1 AU, E so i ~ 10 g cm" 2 (|Weidenschillindl 1977bT ) . Thus, 
at the location of the Earth this yields mj ^ 10 -9 Earth 
masses. 

This shows that even in the most optimistic scenario one 
must start with very small bodies a nd continue by buildin g 
up larger bodies in collisions (e.g., iMorbidelli et al.l I2009T ). 
The origin of the problem here is in the differential shear 
of the disc that limits the l argest possible wavelength of 
gravitational instability (e.g., |Toomrg |l_964) to only ~ 10 8 
cm (see §7. 31 below). 

Now, the mass of the solids inside of a bound fragment 
in the gravitationally fragmenting gas disc is much larger, 
e.g., 
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which for H/R ~ 0.1 yields at least a few Earth masses. This 
is because (a) the parent gravitationally collapsing body is 
a gas disc, not a dust one, and (b), the Keplerian shear at 
100 AU is 10 3 times weaker than it is at 1 AU. 



7.3 Kepler shear and turbulence: preventing 
collapse 

Kepler shear has another undesirable effect on the potential 
gravitational collapse of solid material inside the protoplan- 
etary disc. To become gravitationally unstable, the density 
of the layer must exceed the tidal density, which is of the 
order of 10 -7 g cm -3 at 1 AU. Thus the thickness of the 
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layer must be less than -H" max = £ so i/pt ~ 10 8 cm. In equi- 
librium, the thickness of the layer of solid particles depends 
on their random, e.g., turbulent, velocities. These must be 
less than vt ^ (-Hm ax /-R)«K ~ 10~ 5 v k = 30 cm s _1 . This is 
extremely small compared with the expected sound speed of 
~ km s _1 . Such small turbulent velocities appear very un- 
likely give n the instabilities in th e turbulent Ekman bound- 
ary layer (I Weidenschillind 1 19801 ) excited by the differential 
rotation of solids and the gas. Thus building up solid bod- 
ies as small as 10 -9 Mq by gravitational collapse in the CA 
model is not straightforward. 

Now contrast this picture to the gravitational collapse 
of the dust sphere settling inside the first core, a process 
analogous to the dust layer gravitational collapse. As shown 
in paper I, gas pressure effects become negligible when 
the grain sphere contracts to the "grain cluster" density 
Pgc = ptcfg ■ Self-gravity of the grain cluster at this point 
exceeds gas pressure gradients. The radial size of the grain 
cluster is R gc = R{ c f g , which is typically about 0.1 AU, 
and the mass is f g M{ c . The turbulent velocity required to 
keep the grain sphere from collapse is 
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where c s = y/kT/fi ~ ■\jGMi c JRf c is the sound speed of 
the first core. As temperature of the first core is of order 10 3 
K, we find 
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This velocity is roughly 30% of the gas sound speed in the 
first core, and is about 2000 times higher than the minimum 
turbulent velocity found for the dust layer in the CA model. 
It seems very doubtful that the first core would support such 
a strong turbulence. In contrast to the differentially rotating 
Keplerian disc, the first core may be expected to be in a solid 
body rotation with angular frequency characteristic of the 
disc at R p ~ 100AU, e.g., 3 orders of magnitude slower than 
at 1 AU. 

These estimates are supported by our numerical results: 
even the highest turbulent parameter run, ad = 0.1 (LOal), 
yielded a rather massive solid core (cf. Table 1). 



7.4 Collisional buildup of planetesimals in CA 
model 

As the direct gravitational collapse to form (even such small 
bodies as) planetesimals is complicated by turbulence in 
the CA model, alternative ideas to form them via non- 
gravitational growt4i_jnficliamsms liaye been thoroughly ex- 
plo red (e.g.. iDominik fc Tielenslll997l) ; for a recent review 
see iBlum fc Wu rm (2008J). One could hope that the solids 
in the disc continue to grow via collisional agglomeration 
as smaller dust is believed to do. However, there are two 
principal problems with this hit-and-stick idea also. Firstly, 
due to drag from the gas, the grains of different sizes are 
calculated to move at different speeds; ~ metre size objects 
shou ld be colliding at velo cities of ~ ten metres per sec- 
ond (Weidcnschilling 1977a). Common sense and lab experi- 
ments ( Blum fc WurmlbOOa ) suggest that growth, i.e., stick- 
ing of the collisions partners, is rather unlikely at these high 



speeds. The second problem is frequently referred to as the 
metre-size barrier. In a non-turbulent protoplanetary disc, 
the ~ metre size objects must spiral radially inward, into 
the parent star, on time scales of order a fe w hundred years, 
well b efore they could grow to larger sizes (|Weidenschilling| 
Il977ah . 



7.5 Planetesimals are not needed to build 

terrestrial planet cores inside giant embryos 

We find that neither the problem with too-energetic colli- 
sions for the hit-and-stick growth nor the radial drift prob- 
lem appear to exist if solid bodies grew inside the protective 
environment of the giant embryo. The radial drift problem 
is most obviously absent: solids migrate with the embryo 
rather than individually. 

Furthermore, we shall now show that solids that con- 
tribute the most to the terrestrial planet cores inside the 
giant embryos are not planetesimals, and not even boulders, 
but are meager pebbles. The collisions of such small bodies 
are gentle. This resolved the issue with too fast collisions 
that tend to fragment dust aggregates rather than bind col- 
lision partners together in the CA model. 

As we found in paper I, while grains are microscopic, 
turbulence keeps them from sedimenting down very effi- 
ciently The equilibrium reached between sedimentation and 
turbulent mixing with the turbulent diffusion coefficient 
ad sets the radial scale of the dust distribution Hd ~ 
Rtc(adPtcRtc\/ PaO 2 ) 1 ' 2 (see paper I). We picked the "large" 
grain size limit here, when a> A~4 cm (10 _9 /pf c ), the 
mean free path in the gas. Requiring now that the equilib- 
rium scale-height Hd be smaller than the grain cluster size 
-Rgc - at which point the grain cluster could collapse due to 
self-gravity - we infer that grains must be larger than 
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where q_3 = ay/lO -3 . 

The estimate above predicts that the grains should grow 
to about 1 m size in order for gravitational instability to 
overwhelm turbulent mixing. However, we must keep in 
mind that the estimate is based on a number of approxi- 
mations; e.g., estimating the properties of the "grain clus- 
ter" based on the maximum gas pressure gradient possible 
(paper I) . Numerical experiments in the present paper have 
shown that the "runaway phase" of core accretion, which 
occurs when the grain cluster starts to collapse, takes place 
when the grains reach dimensions of 10-20 cm (cf. figures in 
this paper), rather than 80 cm. 

We can now estimate the sedimentation velocity of 
grains of this characteristic size: 
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Since Rf c /c s « l/^/Gpt c , we can rewrite this as a function of 
the gas density of the first core at the moment of the grain 
cluster formation: 
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This estimate is consistent with our numerical results - the 
grain sedimentation velocity is below 10 m/s in the inner- 
most part of the first core identifiable as the "grain cluster" . 
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One caveat here is that we have so far assumed a 
non-rotating giant embryo. However, comparing the stop- 
ping time of a grain defined as u a /\du a /dt\, one finds that 
only grains larger than 10 3 cm or more would be in the 
"perturbed Kepler orbit" regime (|Weidenschillinelll977al ) in 
which the particle moves essentially on its Keplerian orbit 
experiencing little drag from the gas. The particle's angular 
momentum suspends it some distance away from the cloud 
centre, and it spirals in slowly due to gas drag forces. This 
implies that small particles dominating the grain sedimenta- 
tion process are in the opposite regime, e.g., approximately 
following the rotation of the gas. As long as the embryo is not 
completely supported by rotation, the local rotation speed is 
lower than the circular orbit value. Thus small grains cannot 
be supported by rotation against gravitational collapse, and 
our spherically symmetric treatment should be reasonably 
representative of slowly rotating giant embryos. 

Having said this, we note a sort of reincarnation of 
the "metre-size barrier" problem for dust in rotating giant 
planet embryos - even larger than 1 metre dust grains are 
going to be lost into the centre of the embryo - but that is 
not a problem at all. The grains sinking into the centre of 
the cloud are not lost into the star but presumably join the 
solid core there, therefore they are not "lost" in any sense. 



7.6 Summary on direct gravitational collapse 

Summarising these points, it appears that a direct build up 
of massive terrestrial-like cores is possible inside the first 
cores (giant planet embryos) due to a favourable (spheri- 
cal) geometry and low shear. One interesting point of the 
present picture is that the building blocks of the solid cores 
of terrestrial planets are grains of ~ 10 cm size rather than 
planetesimals of km or more size. While numerical 3D sim- 
ulations, well beyond the scope of our paper, are needed to 
model the amount of mass put into bodies of different sizes, 
there is no clear theoretical need to put much mass into the 
planetesimal population inside the embryo at any time. This 
is in contrast to the core accretion model where terrestrial 
planets are made of planetesimals. We hence expect (pend- 
ing future confirmation by numerical simulations) that the 
relative amount of mass in asteroids and comet-like bodies 
compared to that in terrestrial cores is much smaller in our 
model than traditionally assumed. 

Finally, there is no loss of solids problem as in the disc 
geometry - any solids lost into the centre of the giant embryo 
actually contribute to the growth of the solid core there. 



8 DISCUSSION 

8.1 Summary of numerical results 

We have argued that clumps formed by the GI in proto- 
planetary discs a t R ~ 100 A U should be similar to that 
of the first cores (|Larsonlll969r ). Using a two-fluid radiative 
hydrodynamics code, we found that the dust inside the gas 
clumps may grow and settle into a core if the gas remains 
cooler than the grain vaporisation temperature. The high-Z 
elements core may become massive (~ 10 M$) if feedback 
processes are not too severe. Our r esults are more restrictiv e 
than analytical estimates of lBossi (1 19981 ) mBqss et al.l (|2002l ). 



On the other hand, the maximum mas ses of the sol i d core s 
that we find are larger t han those of iHelled et al.l (|2008f ) ; 
iHelled fc Schubertl (J200Sft . 

As already indicated in the Introduction, the main dif- 
ference with the latter authors appears different initial con- 
ditions for the gaseous clumps. For e xample, all of the 
M > 5Mj mass clumps in Table 1 of IHelled fc Schubertl 
(|2008r ) are unable to form any solid cores because they are 
already hotter than T — 1400 K at birth. In contrast, the 
initial temperature of a M{ c — 5Mj clump in our model 
is only about 60 K (cf. equation [4]). These clumps do cool 
initially very rapidly, increasing in temperature to hundreds 
of K, but their increased density allows for an accelerated 
grain growth too (see fig. 1-3 in paper I). We also found 
that our initial conditions are stable to convection, and thus 
convective grain mixing does not significantly impede grain 
sedimentation. Convection becomes important once a mas- 
sive and radiatively bright solid core forms, but by that time 
grains may be large and the central region dominated (by 
density) by the grains, thus the effects of convection remai n 
subdued compared with the results of IHelled et al.l (|2008l ): 
IHelled fc Schubertl J2008I ). 



8.2 Astrophysical implications 

We discuss observational implications of these results in a 
follow-up paper (paper III) where we take into account ad- 
ditional processes not discussed here: radial migration of the 
giant planet embryos and their Tidal and Irradiative Evapo- 
ration (TIE) when the embryos come too close to the parent 
star. Therefore here we only briefly mention the implications 
directly following from the results presented here and in pa- 
per I. 

We believe that grain sedimentation inside the giant 
planet embryos at tens to a few hundred AU distances from 
the parent star is a viable process by which to form terres- 
trial mass solid cores. We argued in !j7|that a direct gravita- 
tional collapse of ~ 10 cm-sized solids in the central region of 
the gas clump can make massive solid cores directly without 
the need to go through objects of intermediate sizes (such as 
planetesimals and Moon-sized objects). Neither turbulence 
nor the high speed collisions (absent in TIE in the stages be- 
fore the grain cluster collapse) appear to be serious obstacles 
to the gravitational collapse of the solid component. 

We also argued that massive gaseous atmospheres 
around the solid cores should start building up once the 
latter reach a critical mass which is in the range of a few 
to a few tens Earth masses. The arguments here are similar 
to the cl assical core accretion picture for b uilding the giant 
planets l|Mizuno|[l980l ; iPollack et alJll996h except that the 
gas is accreted not from the disc but from the surrounding 
gaseous envelope of the first core. This potentially can make 
a metal rich (as the interior of the first core is such) giant 
planet with a solid core. Removal of the outer metal poor 
gaseous envelope, discussed in the follow-up paper, could 
then produce classical giant planets such as Jupiter. 

In addition to the tidal disruption and irradiative evap- 
oration of the envelope by the parent star, gaseous clumps 
may self-unbind themselves due to an excessive heat re- 
lease by the solid cores within them if envelope opacity is 
high. These effects may be even more important for rotat- 
ing first cores (rotation is expected for first cores formed 
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in a disc, e.g.. iLevinl 120071 ). Indeed, the amount of energy 
needed to unbind the core is then reduced. Furthermore, 
following the suggestion of iBoss et al.l (|2002l ). the gaseous 
envelope of the first cores may be removed by a strong ion- 
ising radiation of a passing OB star. This is clearly more 
realistic for f irst cores as these are initially very fluffy ob- 
jects (see also lWhitworth fe Zinneckerll2004T ). but more work 
is needed here. We plan to address the issue in our future 
work. 

The processes discussed here may create terrestrial mass 
solid planets at large distances from the parent star. In ad- 
dition, the escape velocity form first cores is less than a few 
km s _1 . Collisions of such rotating clumps with each other 
may perhaps unbind a fraction of the solid material there, 
contribu ting to the population of solids at these large radii 
(see also lClarke fc Lodatdl2009t ). 

Finally, we note that the proposed route to formation 
of terrestrial and giant planets is somewhat a of a hybrid 
between the two current competing models, e.g., the grav- 
itational instability and the core accretion models. In TIE 
the planet formation process begins with gravitational insta- 
bility (gas clump formation, grain sedimentation and core 
growth) but continues with as in core accretion (accretion 
of gas onto the solid core) . 



8.3 Shortcomings of the calculations 

Despite using a relatively sophisticated two-fluid radiation 
hydrodynamics approach, our treatment here is still quite 
far from being sufficient. In particular, hydrogen molecules 
dissociation and hydrogen ionisation should be taken into 
account in the future to allow higher resolution studies of 
the innermost region near the solid core. Introducing sev- 
eral different species of grains with different vaporisation 
temperature is desirable, as is allowing the grain size to 
vary from bin to bin, instead of assuming one size as we 
have done here. Turbulent and convective mixing for grains 
may be improved. Furthermore, we have neglected effects 
of rotation of the first core here. Finally, the dust opacity, 
independent and fixed in our calculation, should of course 
evolve as dust grows and sediments, but it is hard to see how 
a self-consistent and yet not overly complicated treatment 
of these effects could be done at the present. 

We expect these effects to modify the results of our 
dust sedimentation calculations. Nevertheless, we feel that 
our exploratory calculations are interesting in showing the 
possible range of behaviour in the given problem. 



8.4 Connection to previous ideas and work 

As spelled out in the Introduction , our work is an ex tention 
and an update of the [Boss] (| 19981 ) ; IBoss et all (|2002T ) ideas. 
We find that giant planet embryos born at R ~ 100 AU 
should be similar in their initial properties to the first cores. 
These cooler and fluffier gas configurations present a perfect 
birthplace for massive solid cores. In paper III we argue 
that giant planet embryos migrate inward and get tidally 
disrupted there. This hypothesis is similar to the ideas of 
Boss and co-authors, although the exact protoplanetary disc 
setting is different. 

There is also some similarity between our ideas and that 



of iRice et all (|2004) ; IClarke fc Lodatol (|2009l ), who have re- 
cently argued that some fraction of planetesimals may form 
inside gas pressure maxima provided by spiral arms in self- 
gravitating protoplanetary discs at R ^ 30 AU. These spiral 
arms are stable for about a rotation, during which time the 
density enhancement in the arms, compared to the mean 
density in the disc, may be significant if the cooling time 
in the disc is short enough. The location of the marginally 
stable non-fragmenting region of the disc is somewhat in- 
side the fragmenting r egion that we discuss (see Fig. 1 in 
IClarke fc Lodatol 120091) . Our study is thus an extension of 
these ideas on the larger disc scales, where spiral arms 
do fragment, and where the density enhancement is much 
greater than unity. 

Furthermore, we note that, having submitted papers I 
and III, and having almost finished the final draft of the 
present paper we discovered a set of even more closely re- 
lated ideas on formation of terrestrial pl anet cores inside g i- 
ant planet e mbryos in a recen t pape r bv lBolev et alj (|2010f l . 
In addition, iBolev fc Durisenl <|2010f ) present numerical sim- 
ulations of fragmenting gaseous disc which includes a popu- 
lation of dust grains and also argue for grain sedimentation 
inside the giant planet embryos. While we have so far con- 
centra ted on the internal evol ution of individual giant em- 
bryos, IBolev fc Durisenl (|2010j) paid more attention to the 
interaction of the embryos with their parent gas disc. These 
completely independent results can be cited as a direct sup- 
port of our ideas. 



9 CONCLUSIONS 

We argued that gas clumps spawned by gravitational insta- 
bility in massive young proto-planetary discs may be excel- 
lent sites for grai n growth an d sedimentat i on, in agreement 
with the ideas of lBosd [|l99Sh ; IBoss et alj (|2002f ). This may 
allow a direct gravitational collapse of the grain component 
of the central part of the giant planet embryo. This process 
spawns terrestrial mass solid cores without going through 
km-sized bodies ($7J. It was also suggested that, similarly to 
the core accretion scenario of giant planet formation, there is 
a critical solid core mass at which the gas atmosphere around 
the solid core collapses to much higher densities. We suggest 
that this process may create metal rich giant planets with 
solid cores inside the massive gas clumps. In a companion 
Letter we show that these clumps should migrate inwards 
where they are disrupted by tidal and irradiative evapora- 
tion, releasing the inner planet. We believe this new hybrid 
scenario for planet formation is promising and should be 
explored further. 
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Table 1. Parameters of the simulations (see £13.61 for detail). ID: id of the simulations. The mass of the first core, M{ c , is in Jupiter's 
masses; fco is the opacity coefficient in the opacity law given by equation [T] a is the power law index in that law; f g is the grain mass 
fraction of the core; t en( j, in IO 3 years, is the duration of the simulation; t acc is the time at which the runaway accretion on the high-Z 
core sets in, in units of I0 3 yrs; M core is the core mass at t = t ent i; M c i and M c 2 are same but at times t = 2 X I0 3 years and t = 10 4 
years, respectively. 



ID 


M tc 


fco 


"d 


u 


''end 


tacc (IO 3 yrs) 


ivl COYC 


M cl 


M c2 


M0a3 


10 


1 


IO" 3 


0.02 


31.5 


9.1 


21.6 


~0 


0.7 


Mlo3 


10 


0.1 


IO" 3 


0.02 


24 


5.0 


5.2 


~0 


4.0 


M2a3 


10 


0.01 


io- 3 


0.02 


10 


1.9 


1.0 


0.2 


1.0 


M3a3 


10 


0.001 


io- 3 


0.02 


1.5 


- 


IO" 6 


- 


- 


M2o2 


10 


0.01 


IO" 2 


0.02 


7.2 


3.0 


0.08 


IO" 6 


0.08 


M2q4 


10 


0.01 


10" 4 


0.02 


6.4 


1.8 


0.7 


0.12 


0.7 


H0o3 


20 


1 


IO" 3 


0.02 


30. 


12.5 


12.2 


3 x 10" 5 


0.0015 


H1q3 


20 


0.1 


IO" 3 


0.02 


10.5 


5.2 


1.5 


2.8 x 10" 5 


1.5 


H2a3 


20 


0.01 


IO" 3 


0.02 


3. 


- 


IO" 6 


IO" 6 


10~ 6 


H3a3 


20 


0.001 


IO" 3 


0.02 


3. 


- 


~ 


~0 


~0 


L0a3 


5 


1 


IO" 3 


0.02 


39.7 


~9.0 


10.8 


2 x IO" 5 


0.15 


Lla3 


5 


0.1 


IO" 3 


0.02 


20. 


6.5 


10.8 


1.5 x 10" 5 


3.4 


L2a3 


5 


0.01 


IO" 3 


0.02 


32 


2.4 


6.2 


6 x IO" 5 


6.1 


L3a3 


5 


0.001 


IO" 3 


0.02 


2.2 


0.7 


1.2 


1.2 


1.2 


LOcvl 


5 


1 


0.1 


0.02 


135 


NA 


0.075 


2 x IO" 5 


10" 4 


L0a2 


5 


1 


0.01 


0.02 


46 


NA 


0.005 


2 x IO" 5 


3 x 10" 4 


Lla4 


5 


0.1 


10" 4 


0.02 


22. 


5.2 


10.8 


5 x IO" 6 


4.4 


L2a2 


5 


0.01 


0.01 


0.02 


14.4 


4.2 


3.1 


2.2 x 10" 5 


2.9 
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